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The Cu40Cl6daca4 system composed of isolated Cu"^^ S=l/2 tetrahedra with antiferromagnetic 
exchange should exhibit properties of a frustrated quantum spin system. Ab initio density functional 
theory calculations for electronic structure and molecular dynamics computations suggest a com- 
plex interplay between magnetic exchange, electron delocalization and molecular vibrations. Yet, 
extensive experimental characterization of Cu40Cl6daca4 by means of synchrotron x-ray diffraction, 
magnetization, specific heat and inelastic neutron scattering reveal that properties of the real mate- 
rial can be only partly explained by proposed theoretical models as the low temperature properties 
seem to be governed by a manifold of molecular configurations coexisting at high temperatures. 

PACS numbers: 75.50.Xx, 36.40.Cg, 63.22.-m 



I. INTRODUCTION 

Strongly frustrated quantum s pin syst ems show, fol- 
lowing numerous theoretical worlpj^isiiisj^ highly degen- 
erate ground states. When the degeneracy is lifted a 
number of exotic ground states may emerge. However, 
in spite of large experimental efforts such examples are 
rare and are hard to probe. Any perturbation tends to 
push the real system into a conventional classical ground 
state. A good recent example is the Cu2Te205X2 (X— CI, 
Br) system. Being a ssumed to be built of weakly cou- 
pled tetrahedral unit^S^^'^'^Iinili] jt appeared to have an 
incommensurate 3D long-range ordered ground stat^^^. 
Our aim is to investigate the properties of tetrahedra 
which are the building blocks of geometrically frustrated 
quantum spin lattices (e.g. kagome, pyrochlorc). Note 
that the ground state of an isolated tetrahedron with 
antiferromagnetic Heisenberg exchange is doubly degen- 
erate in the cas e of h igh T^; symmetry. It is formed 
by two entangledSnil singlets; one being a product of 

two dimers = ^(ITni>+ UiTT)- ITiiT)- liTTi)) 

and the second composed by two-site triplet products, 

= 712(2 |Tin) + 2 litit)- ITUi)- ITUT>- lUTT 

)^ UTTi))- Three triplets and one quintet constitute 
the excited states. We choose in the present work a 
simple frustrated quantum spin candidate among zero- 
dimensional systems. 

We focus our study on the CU4OCI6L4 system. Its 
main magnetic building block is an almost ideal isolated 
tetrahedral cluster of four Cu^"*" ions (see Fig. [1]). The 
magnetic properties of several representatives of this fam- 
ily could not be quantitativel y accoun ted for by a simple 
Heisenberg isotropic exchang^i^ESIiZl A number of more 



complex models have been proposed. In early stages 
Lines et al^ suggested the existence of orbital degener- 
acy of the Cu^^ ion and a large antisymmetric contribu- 
tion to the spin Hamiltonian. However, angular overlap 
model calculations'^ indicated that the unpaired electron 
of the Cu^"*" ion occupies the dz^ nondegenerate orbital 
with no orbital moment, which is around 960 meV above 
other d-orbitals. 

Henceforth the invoked models follow two routes, the first 
one exploits the complexity of the spin exchange and 
the second one puts forward the spin-vibrational inter- 
actions. The first approach is based on the generalized 
spin Hamiltonian 

Hex -^tso ^" Has ^" -^an ^ ^ Jij^i^j ^" 

^G„-(S,xS,)+^ J^SfS^" (1) 

ij ij a—x,y,z 

where Hiso is the isotropic Heisenberg-Dirac-van Vleck 
term, and Has and Han are the antisymmetric 
and a nisotrop ic terms respectively. A number of 
author^i^ESlsni explained magnetic susceptibility, elec- 
tron paramagnetic resonance (EPR) and inelastic neu- 
tron scattering (INS) observations by considering large 
antisymmetric contributions (up to 20 % of the isotropic 
exchange constant). Dickinson et al^ suggested the 
presence of ferromagnetic intercluster exchange in addi- 
tion to the antiferromagnetic intracluster one. 

The second approach is based on the spin-vibrational 
Hamiltonian 

H^H,,b{q) + H,,iS,q) (2) 

where Hyii,{q) describes the harmonic vibration of the 
nuclei around the equilibrium position and Hex{S,q) 
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depends not only on the spin S but also on the nuclear 
displacements q of the Cu^"*" ions. JoneJ^^ developed a 
simplified model of dynamic and fluxional Jahn- Teller 
distortions of CU4OCI6L4, assuming a linear relation 
between the coup ling constants and atomic shifts. 
Polinger et al^^^ extended this approach concluding 
that the pseudo Jahn-Teller effect is the most probable 
reason for the anomalous magnetic properties of the 
system. 

We collected substantial experimental data on 
CU4OCI6L4 with L=diallylcyanamide=N=C-N-(CH2- 
CH=CH2)2 (Cu40Cl6daca4) and performed ah initio 
density functional theory calculations for the electronic 
structure as well as ah initio molecular dynamics. Our 
findings suggest that indeed the vibrational degrees of 
freedom are very important. However, the dominant 
ingredient is the presence of a manifold of molecular 
configurations coexisting at high temperatures. There- 
fore both, the generalized spin and spin- vibrational 
Hamiltonians remain insufficient to explain the mea- 
sured magnetic properties of CU4OCI6L4. 



II. EXPERIMENTAL DETAILS 



A has been used for comparison with the FOCUS results 
on the protonated sample. Empty sample holders have 
always been measured and subtracted. 



III. CALCULATION METHOD 

A. Electronic structure 

In order to understand the electronic and magnetic 
properties in this compound, first principles local spin 
density aproximation calculations were done within the 
framework of the Density Functional Theory (DFT). For 
the description of the exchange and correlation energy, 
the Generalized Gradient Approximation (GGA) was for- 
mulated by the Perdew-Burke-Ernzerhof (PBE) density 
functionapS. 

In the spin-polarized version of the GGA approximation 
the exchange-correlation energy E^^c is a functional of the 
local electron spin densities n| and n^ and their gradients 



The Cu40Cl6daca4 single crystals (Nr. 1 and Nr. 2) 
and the polycrystalline sample have been prepared by 
reaction between CuCl and diallylcyanamide in the pres- 
ence of 1J2II. The crystal structure has been examined by 
X-ray synchrotron powder and single crystal diffraction 
performed at the BMOlA Swiss-Norvegian Beamline 
of the European Synchrotron Radiation Facility ESRF 
Grenoble, France at a wavelength of 0.7272 A in the 16 
K - 340 K temperature range. 

The bulk properties have been studied using a Quantum 
Design PPMS commercial system. The dc- and ac- 
magnetic susceptibility of polycrystalline and single 
crystal materials have been measured in the 1.8 K - 
300 K temperature range and in applied magnetic fields 
up to 9 T. The measurements have been performed 
mostly on cooling with the cooling speed of 0.5 deg/min. 
The specific heat has been collected in 0.35 K - 10 K 
temperature range and in applied magnetic fields up to 
14.5 T, and the magnetization - in the 1.8 K - 20 K 
temperature range up to 14.5 T. 

The inelastic neutron scattering experiments have 
been performed on a 2 g polycrystalline protonated 
Cu40Cl6daca4 sample on the neutron time-of-flight 
spectrometer FOCUS at the spallation neutron source 
SINQ, Paul Scherrer Institute, Switzerland and on a 
0.5 g polycrystalline deuterated sample on IN4 at the 
Institut Laue Langevin ILL Grenoble, France. Four 
setups of FOCUS have been exploited {\i= 5.45 A, 
4.87 A, 3.42 A and 2.26 A) to access the - 13 meV 
energy-transfer range with various energy resolutions. 
On IN4 high-energy transfers up to 40 meV have been 
probed with Ai=1.32 A, while a setup with Ai=2.64 



E^^^[n^,n^]= f{nT,ni,VnT,Vni). (3) 



For calculating the eigenvalues and eigenvectors in the 
DFT scheme the Full Potential Linearized Augmented 
Plane Waves method (FP-LAPW), as implemented in 
the Wien2k cod^^, has been used. In this method the 
unit cell is divided in two regions, the first one is com- 
posed by non-overlapping spheres centered at the atomic 
sites and the second is the region between the atomic 
spheres or interstitial region. In the first region a spheri- 
cal potential is assumed and therefore the wave functions 
are expanded in terms of spherical harmonics. In the in- 
terstitial region a constant potential is assumed. This 
allows to expand the wave functions in that region in 
terms of plane waves. 

In our calculation the valence wave functions have been 
expanded up to ^=10. Sphere radii Kmt of 2.1 a.u. for the 
Cu atom, and 1.6, 1.4 , 1.2, 0.9 and 0.84 a.u. for CI, O, 
N, C and H atoms respectively have been set. Within the 
Cu muffin-tin radii chosen we excluded the 4s states (the 
average radii for the wave functions belonging to these 
states is 2.7 a.u.), therefore these states lie in the inter- 
stitial region. The Cu-3c? and 3p states enter in the cal- 
culation as the valence and semicore states respectively. 
For the Cu 3c? states the LAPW basis with augmented 
plane wave plus local orbital [APW-I-LO] method have 
been used. The plane-wave expansion with Hmt^max 
equal to 3.76 has been set. This is a typical value when 
C-H bonds are present. For the integration in the irre- 
ducible wedge, 48 k-points have been used. 
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B. Molecular dynamics IV. RESULTS 



The dynamics of Cii-daca has been investigated using 
the DFT code VASP^^^s ^i^i^h is well adapted to 
a system as large as Cu40Cl6daca4 (174 atoms per 
unit cell). VASP uses a plane wave basis set and the 
Projector Augmented Wave (PAW) pseudo-potentials to 
describe the core electrons of the atoms. The GGA-PBE 
functional was usecP^ and all electronic calculations 
were performed at the T point in reciprocal space. Since 
phonons and molecular vibrations are of interest in 
this paper, lattice dynamics calculations, involving the 
determination of the dynamical matrix, were attempted. 
In view of the number of atoms in the unit cell, the 
direct method^^ is the only practically feasible approach 
in which symmetry-inequivalent atoms are displaced one 
by one from equilibrium and the Hessian is constructed 
from the Hellmann-Feynman forces from the series of 
electronic calculations. 

In the case of Cu40Cl6daca4, the direct method always 
resulted in negative frequencies. This method invokes 
a significant perturbation to the equilibrium configu- 
ration and therefore the electronic structure. Indeed, 
electronic and structural correlations between the core 
and the ligands lead to an electronic structure which 
is very sensitive to the precise molecular configuration. 
Such sensitivity is an indication of the fragility of this 
molecular system as proved by the physical meaning of 
negative frequencies in this case. 

In this situation we resorted to ab initio molecular 
dynamics, using VASP, in which case the (partial) vibra- 
tional density of states ((p)VDOS) can be obtained from 
the velocity auto-correlation function and the scattering 
function is obtained from the van Hove correlation 
functions. Both of these correlation functions are deter- 
mined directly from the ah initio molecular d yna mics 
(AIMD) trajectories using the nMoldyn prograrrPSl. The 
temperature dependence of the dynamics, which may 
reveal anharmonicities, is investigated via the average 
kinetic energy of the atoms. The drawback of the 
molecular dynamics approach compared to a lattice 
dynamics calculation, is that the phase relation of the 
atomic displacements is lost and normal modes and their 
frequencies cannot be determined. One solution is to 
investigate the time-dependent fluctuations of geometric 
quantities associated with groups of atoms. In the case 
of the magnetic properties of Cu40Cl6daca4 the Cu 
tetrahedron is of interest. 

AIMD simulations have therefore been performed at 100 
K and 250 K for a protonated crystal. Production simu- 
lations at constant volume and energy (NVE ensemble) 
cover 10 ps, that is KiW^ simulation steps. The 100 K 
simulation was then re-equilibrated for a deuterated crys- 
tal and the NVE production run in this case lasted 27 ps 
in order to give better access to the low frequency modes. 



A. Crystal structure 

The Cu40Cl6daca4 structure is tetragonaP^ (space 
group P42ic) in the whole 16 K - 340 K studied temper- 
ature interval. There are two Cu40Cl6daca4 molecules 
per unit cell and they are slightly twisted around the z 
axis (Fig.g. The oxygen is in the center of the molecule, 
surrounded by a tetrahedron of Cu^"*" ions with two differ- 
ent Cu-Cu distances, equal to 3.0790(2) A and 3.1185(2) 
A at 80 K. The tetrahedron is enclosed by a distorted 
octahedron of chlorine ions with four Cli in the square 
base and two CI2 in the apexes. Each copper has trigonal 
bipyramidal coordination: three Cl^ ions form the trian- 
gular base of the bipyramid, while the central ion 
and a terminal ligand L are located in apical positions. 
The exchange paths Cu-O-Cu are identical, while the two 
paths through chlorine ions, Cu-Cli-Cu and Cu-Cl2-Cu, 
slightly differ (see Table |l]). Such geometry implies rather 
high symmetry of the cluster and of the spin exchange 
Hamiltonian, only marginaly deviating from T^. 
There is an anomaly in the lattice constants in the 230 

FIG. 1: (Color online)The ab projection of the Cu40Cl6daca4 
unit cell. The Cu atoms are represented by red circles, CI, N 
and C - by green, blue and grey circles, respectively. H atoms 
are omitted for clarity, O atoms are located behind C12 and, 
therefore, are not visible. 




K - 280 K temperature range notifying an isostructural 
order-disorder crossover^. Apparently the molecule at- 
tains several configurations with small energy difference 
between them. Above Tc=282 K these configurations 
coexist, while below Tc molecules freeze in one of them. 
The configurations differ mainly by the position of the 
branched diallylcyanamide ligands, which fold. This 
modifies the distances between the molecules (Table |l]) 
and the unit cell volume shrinks by 7 % at 80 K compared 
to the volume at room temperature. The CU4OCI6N4 
core changes very little. The only noticeable difference 
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TABLE I: Comparison of the characteristic angles [deg], 
intra- {di„tra [A]) and inter- (dinter) molecule distances in 
Cu40Cl6daca4 at 80 K and 300 K. Mean-square displacement 
amplitudes of Cu (MSDcu, [A^]) and other atoms (MSD;ig, 
[A^]) in the specified bond directions of the CuClaON trig- 
onal bipyramid calculated from the anisotropic displacement 
parameters refined at 80 K. 



-J 


oU IV 


OUU XV 


Angles 


80 K 


300 K 


Cu-Cu 


3.0790(2) 


3.0933(4) 


'^CuOCu 


108.106(6) 


108.9(2) 


Cu-Cu 


3.1183(2) 


3.1101(5) 


^CuOCu 


110.158(6) 


109.8(2) 


Cu-O 


1.9016(2) 


1.9011(3) 


^NCuO 


179.01(4) 


178.8(1) 


Cu-Nl 


1.932(1) 


1.935(3) 


^CuCllCu 


80.46(1) 


80.48(3) 


Cu-Cll 


2.3694(4) 


2.3911(9) 


^Cucncu 


80.59(1) 


80.62(3) 


CU-C12 


2.3806(3) 


2.3909(8) 


■^CllCuCll 


117.97(2) 


118.85(4) 


Cu-Cll 


2.4578(3) 


2.4231(9) 




127.73(1) 


120.92(3) 




MSDc« 


MSD,,,, 


dinter 


80 K 


300 K 


Cu-O 


0.0107 


0.0125 


dlOO 


12.60490(3) 12.5652(1) 


Cu-Cll 


0.0143 


0.0130 


dool 


12.73045(5) 13.7685(2) 


CU-C12 


0.0203 


0.0189 


dm 


10.9525(3) 


11.2399(1) 


Cu-Nl 


0.0108 


0.0152 









is the shift of the Cli ions by 0.2487(9) A, while for Cu 
and CI2 the shift is only 0.0462(9) A and 0.0282(9) A, 
respectively. 

Inspection of the interatomic distances and angles pre- 
sented in Table |T] suggests that the Cu-O-Cu magnetic 
exchange path is antiferromagnetic (AF), while the Cu- 
Cl-Cu one is ferromagnetic (F). Therefore the mobility 
of the Cli ions is essential in determining the magnetic 
exchange and this is investigated with the help of AIMD 
simulations. 

One more appreciable detail could be extracted from the 
single crystal diffraction data. The difference Fourier 
map through the N-O-N plane of the molecule (Fig. [2]) 
shows additional electron density near Cu atoms at 80 K. 
The tiny peaks (0.38 and 0.31 e/A'^) are displaced along 
z by 0.84(5) A and 0.73(10) A from Cu towards O and 
N (Figure [2]) , respectively. At 300 K these peaks are less 
pronounced (0.24 e/A'^). These difference electron den- 
sity peaks near Cu and O, N, C positions evidence the 
delocalization of e-density in the 0-Cu-N=C bonds. No 
smearing of the electron density near copper atoms or an 
increase of their anisotropic displacement parameters is 
observed; the mean-square displacement amplitudes are 
very small (see Table . This implies the absence of a 
static and/or dynamic Jahn- Teller effect at low temper- 
atures. 



B. Electronic structure and charge density 
distribution 

The electronic Density of States (DOS) has been cal- 
culated for the crystal structures determined at tempera- 
tures, 80 K and 340 K. Since both structures have almost 
indistinguishable DOS features, we present here the re- 



FIG. 2: Difference Fourier map of Cu40Cl6daca4 at 80 K 
(top) and 300 K (bottom) from x-ray single crystal data. 
6 r J. - ■ rr-^ -r- ■ — — T-n rn 
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suits for the T=80 K structure. In Figure |3]we show the 
partial Cu d, CI p, N p, O p and C p DOS in the range 
of energies between -4 and 4 eV obtained from non-spin 
polarized GGA^ calculations. We observe a very narrow 
and sharp contribution of Cu d and CI, N and O p states 
at the Fermi level well separated from the rest of the va- 
lence and conduction states. These features indicate that 
these states strongly hybridize and that the system is al- 
most zero dimensional, formed by well isolated Cu tetra- 
hedra. The bandwidth at the Fermi level is less than 0.2 





eV (see the inset of Fig. [s] where a blow up of the DOS 
at the Fermi level is shown). From calculations of the 
crystal field splitting we conclude that the Cu d band at 
the Fermi level has (1^2 character (in the local reference 
frame of the CUCI3ON trigonal bipyramid). 

Within the GGA approximatiorpS we obtain that the 
Cu d states are partially occupied and therefore the sys- 
tem is metallic. This is a well-known drawback of DFT 
calculations within the LDA or GGA approach when 
performed for correlated systems. Improved exchange- 
correlation functionals like LDA+U should provide the 
correct insulating behavior of the material. In the case of 
Cu40Cl6daca4, calculations at the level of spin-polarized 
GGA already open a small gap at the Fermi level as 
shown in Fig. |4] In this figure we present the partial 
DOS for Cu d states in both spin up (upper pannel) and 
spin down (lower pannel) channels. The occupation of 
the spin-up and spin-down Cu d states is very similar, 
therefore the resulting localized magnetic d-moment is 
small. Table [iT] presents the magnetic moment at the 
cluster core inside the muffin-tin radii. The unpaired Cu 
electron is strongly delocalized via hybridization mainly 
with oxygen and less strongly with chlorine and nitrogen. 
In Fig. [5] we show the projection on the xy plane of the 
calculated charge density map of the band around the 
Fermi level. We note the good agreement with the mea- 
sured difference Fourier map from the 80 K x-ray diffrac- 
tion experiment (Fig. |2|. 
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3: (Color online) Partial Density of States in the non- 
polarized GGA calculations for Cu40Cl6daca4 at T=80 



FIG. 4: Partial Density of States for majority (upper panel) 
and minority (lower panel) spin contributions from Cu-d 
states of Cu40Cl6daca4 at T=80 K. 
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C. Bulk properties 



The temperature dependence of dc-magnetization 
(M/H) measured for the polycrystalline sample at H=0.1 
T and for a single crystal (Nr. 1) at several values of ap- 
plied magnetic field is presented in Fig. [6] At 300 K 
the (M/H)-T value is 0.375 cm3(molCu)-^ suggesting 



TABLE II: Magnetic moment (M[ms/atom]) within the 
muffin-tin radii (rmt) and its total moment fraction [%]. For 
the interstitial the moment is normalized per formula unit 
(f.u.=CuOo.25Cli.5daca). 



Atom 
Cu 
O 
CU 
C12 
Nl 
interstitial 



2.1 
1.4 
1.6 
1.6 
1.2 



M 
0.58350 
0.44204 
0.08352 
0.09027 
0.02857 
0.14477 



58 
11 
8.3 
4.5 
2.9 
14 



that the system is composed of uncoupled S=l/2 spins. 
The Curie- Weiss temperature estimated from the high- 
temperature susceptibility is 9i=—50 K, which implies 
dominant antiferromagnetic interactions. 
For the polycrystalline sample M/H gradually increases 
with temperature lowering and at 1.8 K reaches 0.02 
cm'^(molCu)^^. This is a very small value, suggesting 
the presence of a tiny amount of free S=l/2 impurities. 
However, the temperature dependence cannot be fitted 
to a Curie-like 1/T contribution, which would grow much 
faster at low temperatures. 



6 



FIG. 5: (Color online) Charge density distribution at T=80 
K. 




For a single crystal, dc-magnctization consists of two ma- 
jor contributions. The first contribution has a bump near 
12 K and can be attributed to isolated tetrahedra with 
antiferromagnetic exchange with J =-1.8 meV resulting 
in a spin-singlet ground state. The second contribution 
increases with temperature decrease and, similarly to the 
polycrystalline case, cannot be fitted to the paramagnetic 
1/T term. The ratio of these contributions is different for 
different single crystals and depends on crystal orienta- 
tion and applied magnetic field. Based on a number of 
dc-magnetization measurements we reckon that the ra- 
tio of these contributions and the sensitivity to magnetic 
field is related not to a specific crystallographic direction 
but to the morphology of crystal. Apparently there is a 
strong growth anisotropy and not a lattice anisotropy. 
Crystals usually grow as platelets with the [110] axis di- 
rected along the shortest edge, and [001] and [1-10] being 
in the basal plane. An example of M/H utmost sensitive 
to the field is presented in inset of Fig. [6] The mea- 
surement is performed with magnetic field applied along 
the shortest edge. With increasing field the bump at low 
temperatures decreases and already at H=0.5 T vanishes. 
For another crystal (Nr. 2) M/H shows no bump along 
the shortest edge even at H=0.2 T (Fig.[7|, though when 
the magnetic field is applied along the middle [001] or 
long [1-10] edges of the crystal, the bump is tractable to 
much higher fields (not shown). 

One more important observation is that the low-field 
(0.05 T) dc-magnetization curves measured for the zero 
field cooled and 0.05 T cooled conditions are different 
(not shown). This implies that cooling in a relatively 
small field is sufficient to change the microscopic state of 
the crystal. 

The field dependence of dc-magnetization is also sample 
and direction dependent (Fig. [Tjinset). For the polycrys- 
talline sample magnetization is concave and does not sat- 



FIG. 6: (Color online) Temperature dependence of de- 
magnetization (M/H) of the single crystal Nr. 1 measured 
at 0.1 T, 0.2 T, 0.5 T, 4 T and of polycrystalline (black solid 
line) Cu40Cl6daca4 measured at 0.1 T (left scale). (M/H)-T 
of crystal and powder measured at 0.1 T (right scale). In- 
set: zoomed 1.8 K - 30 K temperature interval of M/H of the 
crystal. 
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urate at 14 T reaching only 0.27 ms/Cu. For the single 
crystal Nr. 2 magnetization attains almost 0.5 tob/Cu 
for the [110] direction, while for other directions it ap- 
proaches values similar to the polycrystalline sample. 
The most plausible origin of our observations is the pres- 
ence of nonuniform intratetrahedral exchange couplings 
within a single crystal and - even to larger extent - within 
a polycrystalline sample. This diversity most probably 
appears due to the manifold of molecular configurations 
above the order-disorder crossover at Tc=282 K, which 
could influence the microscopic state of the system in two 
ways: i) during the crystal growth, performed at room 
temperature, several non-ground state molecular config- 
urations freeze-in, ii) during the field cooling through 
the crossover region the applied field favors certain non- 
ground state molecular configurations. In these frozen-in 
molecule configurations the Cli ions could have positions 
which strengthen the ferromagnetic exchange path be- 
tween two or more Cu^"*" ions. As these molecular con- 
figurations are minority and do not represent the ground 
state, we shall call them 'defects' for short. 
In agreement with the above considerations it is reason- 
able to assume a coexistence of at least two kinds of clus- 
ters to explain the measured dc-magnetization: i) AF 
ones, with antiferromagnetic (AF) exchange of the order 
of J=-1.8 meV, ii) AF/F ones, with ferromagnetic (F) 
couplings of the similar strength as the AF ones. The sus- 
ceptibility calculated for this case (Fig. |8] red curve) fits 
relatively well the experiment. We can fit our data also 
with the model comprising large anisotropy of exchange 
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FIG. 7: (Color online) Temperature dependence of de- 
magnetization (M/H) of the single crystal Nr. 2 measured 
with magnetic field H=0.2 T applied along three crystal- 
lographic directions: [110] (green circles), [001] (red trian- 
gles) and [1-10] (blue squares), corresponding to short, mid- 
dle and long edges of the crystal. Inset: Field dependence 
of dc-magnetization of single crystal Nr. 2 and powder 
Cu40Cl6daca4 at 1.8 K. 
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FIG. 8: (Color online) Experimental dc-magnetization (M/H) 
of the single crystal Nr. l(black) and the calculated suscep- 
tibility for i) two types of clusters, AF and AF/F with two 
AF Ji=-1.8 meV, four F J2=1.8 meV (red), ii) one type of 
tetrahedral clusters with Jii,o=-l-8 meV, Jz = 1.2 meV (blue). 
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within the single type of clusters (Fig.jsjblue curve), but 
an unreasonably large anisotropy parameter Jz—1.2 nieV 
is required. We therefore give preference to the hypoth- 
esis of the presence of two kinds of clusters. 

Within the proposed scenario it is also easy to un- 
derstand the strong sensitivity of the bulk properties to 
the applied magnetic field which is significantly smaller 
than J. The real part of the ac-susceptibility of poly- 
crystalline sample shows several steps as a function of 
applied field (Fig. |9|. The most pronounced step is be- 
tween 2 T - 4 T, which corresponds to the Zeeman energy 
of 0.24 meV - 0.48 meV. Also the low temperature spe- 
cific heat of polycrystalline and single crystal samples 
contains a broad Schottky-type anomaly centered at 0.4 
K (0.034 meV) in 0-field data (Fig. flO]). This feature 



shifts to higher temperatures with an applied field. The 
magnetic entropy extracted below 10 K after subtraction 
of a constant phonon contribution reaches only 0.7-0.8 
J/KmolCu. This is approximately 1/8 of the Rln2=5.763 
J/KmolCu value expected for the spin degree of freedom, 
so this magnetic contribution corresponds to a minor part 
of the sample. We deduce that one of the triplet E-levels 
of the AF/F clusters is situated at 0.4 K, though collected 
data are insufficient to propose the full energy scheme for 
AF/F clusters. 



For the protonated sample at 2 K two sharp features 
have been detected at 1.8 meV and 2.4 meV on the neu- 
tron energy loss side, they are labeled 1 and 2 in Fig. [Tl] 
The peaks are at least two times broader than the resolu- 
tion function (0.08 meV at 1.8 meV). At higher tempera- 
tures (7 K, 15 K, 35 K) the 1.8 meV feature broadens and 
shifts to lower energies or a 2nd peak at slightly lower en- 
ergy of 1.6 meV appears and increases in intensity, while 
the 1.8 meV peak decreases. The 2.5 meV feature in- 
creases in intensity and an additional feature at 1 meV 
(labeled 3) appears. No extra peaks have been observed 
in the low energy region down to 0.4 meV. 
A number of peaks are observed at higher energy trans- 
fers (not shown in Fig. Ill, of which the most pronounced 



is the one at 4.8 meV. Its width is at least three times 
broader than the resolution function (0.3 meV at 4.8 
meV) and intensity increases from 2 K to 35 K. 
The measurements on the deuterated sample reveal es- 
sentially the same excitation spectrum (Fig. 12 1. The 



1.8 meV peak tends to vanish at 75 K, while the 2.4 meV 
and 4.8 meV features decrease from 2 K to 75 K. From 
the Q(momentum)- and T(temperature)-dependence we 
relate the 1.8 meV peak to a spin singlet-triplet transi- 
tion between the cluster levels, though it is not clear if 
it is a single excitation or several ones. It clearly cor- 
responds to the bump in M/H. The 4.8 meV feature is 
increasing with Q and most probably is related to molec- 
ular vibrations, while for the 2.4 meV and other features 



8 



FIG. 9: (Color online) Real part of ac-susceptibility of poly- 
crystalline Cu40Cl6daca4 measured in an applied magnetic 
field up to 9 T at 1.8 K, 3 K and 6 K. 




the FOCUS and IN4 results are less conclusive. 
The measurement at higher energy transfers showed that 
relatively sharp peaks discussed above reside on a very 
broad feature (up to 20 meV), which increases in intensity 
with temperature and Q (Fig. 13 1. These spectral fea- 



tures and their magnetic or phonon origin are understood 
with the help of AIMD simulations (below) , which repro- 
duce the inelastic scattering response and allow deeper 
insight via partial vibrational densities of states and spe- 
cific two, three or four atom correlation functions. 



E. Molecular dynamics 

AIMD simulations give a direct view of the amplitude 
of atomic displacements. While the light atoms tend to 
have larger displacements than the heavier atoms (CI and 
Cu), this is not the case for O which is trapped in the 
Cu-tetrahedron and the Cl-octahedron. Note that the 
atoms are classical particles in the AIMD simulations so 
that their displacements do not include zero-point mo- 
tion, which can be of the order of 30% for hydrogen 
at room temperature. Calculated mean-square displace- 
ments (MSD) averaged over all atoms of the same type 
are 0.0177 for Cu, 0.0300 for CI, 0.0093 A^ for O 
and 0.0283 A^ for N. This is in good agreement with val- 
ues obtained by projection on specified bond directions 
from the 80 K crystal structure refinement listed in Ta- 
ble H 

The MD simulations are validated by comparing calcu- 
lated and measured scattering functions. We note that 



FIG. 10: (Color online) Log/log plot of specific heat of poly- 
crystalline (top) and single crystal (bottom) Cu40Cl6daca4 in 
the 0.28 K - 10 K temperature range with applied magnetic 
field T - 12 T. 
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the low frequency modes are independent of hydrogen 
or deuterium (H/D) substitution and must therefore be 
delocalized modes involving many atoms. An acoustic 
phonon would show a 2% frequency shift upon deuter- 
ation. Fig |14| shows the coherent scattering function 
Scoh{Q,^) from the 27 ps simulation on the deuter- 
ated system. Spectra at lower Q-values have weak in- 
tensity and are not reliably calculated with our model. 
Low frequency, numerical noise in the vibrational Den- 
sity of States (vDOS) contributes strongly to the coher- 
ent structure factor since this quantity is proportional to 
vDOS/frequency. At high Q there are peaks at 2 meV, 
4.5 meV and 7 meV and a shoulder at «10 meV which 
are consistent with the experimental data (Figs. TI]|l3 1. 
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FIG. 11: (Color online) INS spectra of protoriatcd 
Cu40Cl6daca4 at T=2 K (blue) and 35 K (red) measured on 
FOCUS with Ai=4.85 A. The intensity has been integrated 
over the momentum transfer ranse 0.5 < Q < 2 A-^ 
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Since the AIMD simulations do not include spin degrees 
of freedom, these peaks can only be the signature of 
lattice excitations and, as expected, their intensities in- 
crease with increasing Q. The coherent scattering func- 
tion is only calculated at the reciprocal lattice points of 
the simulation cell so that Q-sampling is limited to the 
Brillouin zone centre when using a single cell for the MD 
simulations. Simulated spectral peaks therefore corre- 
spond to F-point modes. 

Atomic contributions to the signal are obtained from the 
partial vDOS (pvDOS) which are shown in Fig. 15 top. 
Each pvDOS has been normalized to the corresponding 
number of degrees of freedom, 3n, where n is the number 
of a particular species. Fairly well-defined peaks exist at 
the lowest frequencies in the low temperature MD simu- 
lations. The delocalized nature of the modes is shown by 
the fact that all atoms have similar pvDOS's below «5 
meV. This result is consistent with there being negligi- 
ble spectral shift upon deuteration. The trapped oxygen 
atom has a characteristic vibration frequency of the or- 
der of 55 meV and this is the only atom that does not 
contribute to the lowest frequency modes. The analy- 
sis of the chlorine atoms has been refined to separate 
the contributions from apical atoms and those in the 
square base around the copper tetrahedron. There is 
a small but significant softening of the apical-chlorine vi- 
brations compared to those of the square-base chlorines, 
average frequencies and standard deviations being, re- 
spectively, 18.8±13.5 meV and 20.3±14.4 meV. Fig. [15] 
bottom shows examples of the lowest frequency part of 



FIG. 12: (Color online) INS spectra of deuterated 
Cu40Cl6daca4 at T=2 K (top) and 75 K (bottom) measured 
on IN4 with Ai=2.64 A for three Q-values 1.5 (blue), 2.5 (vi- 
olete) and 3.5 (red) A~^, with increasing intensity integrated 
over 0.5 A"^ 
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the vDOS for D, C and Cu for different lengths of sim- 
ulations. The dotted curves show the results from a 10 
ps trajectory while the solid curves are derived from the 
27 ps trajectory. The longer trajectory gives better fre- 
quency resolution which shows that the first peaks in the 
DOS occur at 2 meV. Spectral peaks below this frequency 
must therefore be of magnetic origin. 
Features of the VDOS can be further understood by con- 
sidering n-atom correlations. The frequency spectrum 
of interatomic distances involving Cu atoms is shown in 
Fig(l6l The Cu-Cu stretch is strongest at «12 meV. The 
twisting deformation of the tetrahedron, that is the an- 
gle variation between Cu-Cu vectors on opposite sides 
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FIG. 13: (Color online) INS spectra of deuterated 
Cu40Cl6daca4 at T—2 K and 75 K measured on IN4 with 
Ai = 1.32 A. The intensity has been integrated for two Q- values 
3.5 and 5.5 A'^ over 1 A'^ 
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of the tetrahedron, also has a characteristic frequency of 
wl2 meV. These are the modes that could directly mod- 
ulate the magnetic copper-copper interactions. The oxy- 
gen pathway is clearly not modified by any low frequency 
modes whereas the copper-chlorine distance modulations 
show a significant amplitude from 25 meV down to the 
lowest frequencies. 

Spin restricted AIMD simulations are the only way to 
gain insight into the low frequency modes measured on 
FOCUS and IN4. Simulations show excitations at above 
2 meV and at ^4.5 meV which grow in intensity with Q 
and therefore correspond mainly to molecular vibrations. 
Below these frequencies, the longest simulations (27 ps) 
indicate that there are no lower lying lattice vibrations 



implying that the lower energy features observed in the 
experiment are of magnetic origin. 

FIG. 14: (Color online) Calculated coherent scattering func- 
tion (CSF) for four Q-values ranging, with increasing inten- 
sity, from 2.5 to 5.5 A~^ of deuterated Cu40Cl6daca4. 
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V. SUMMARY 

The Cu40Cl6daca4 system comprises a very interest- 
ing magnetic unit - a Cu^'*' S=l/2 tetrahedron. Due 
to the geometry of the cluster two exchange paths are 
possible: antiferromagnetic Cu-O-Cu via the central 
oxygen atom and ferromagnetic Cu-Cl-Cu via the 
peripheral chlorine atoms. In the ideal system the AF 
exchange with J=-1.8 meV dominates. 
Our DFT calculations show that unpaired electron 
is only partially localized in the Cu d-states, and 
hybridizes strongly with the oxygen atom in the center 
of the tetrahedron and is further delocalized around 
atomic cores within the molecule. Besides the O^^ 
and Cl^ ions, which certainly participate in the Cu-Cu 
superexchange, the N=C bond might play an important 
role of the charge reservoir. We plan to confirm this 
picture by polarised neutron spin density mapping and 
to study the interplay between electron delocalization 
and magnetic interactions in more detail. 
The presence of strong molecular vibrations emerge 
from our INS experiments and spin restricted AIMD 
simulations. Although all atoms participate in the 
vibrations, it is the organic ligands and the apical Cl^ 
ions of the cluster which vibrate most strongly, while the 
O^^ ion remains pinned in the center of the tetrahedron. 
Therefore the spin- vibrati onal H amiltonian suggested by 
JoneJ^ and Polinger et alJ ^^ ^ ^^ I might be the appropriate 
model to describe the ideal Cu40Cl6daca4 system. 
However, the properties of real samples, cannot be 
that easily attributed to this model. Measurements 
of the bulk properties of single crystals reveal strong 
anisotropy which dominates also the bulk properties 
of polycrystalline samples. We infer that only part of 
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the clusters retain the spin singlet ground state, while 
another part has magnetic ground state and the ratio 
between these fractions is sample and thermal history 
dependent. 

The anisotropy originates most probably from defects 
introduced during the sample synthesis as there is only 
small energy difference between several molecular con- 
figurations. The isostructural order-disorder structural 
crossover at Tc=282 K complicates the situation even 
more, as applied magnetic field and thermal history 
may affect which molecule configuration freczcs-in 
determining in this way the magnetic exchange at 
low temperatures. In this sense the high-temperature 
manifold of molecular configurations governs magnetic 
exchange at low temperature. 

It is well possible that molecular vibrations influence the 
low temperature magnetic properties as well, we leave 
this question open for future investigations. We suspect 
that unusual magnetic susceptibility meas ured for other 
representatives of the CU4OCI6L4 systerrP ^ I ^^ I ^^ I might 
have similar origin to Cu40Cl6daca4. It is important 
to note that influence of the structural and vibra- 
tional degrees of freedoms on magnetic properties of 
molecular magnets is recognized in a number of other 
systems, for example in {V15} comple3p^,{A^j4Moi2} 



spin tetrahedrorPsl^ {CU3} spin triangl^^. 
Finally, the complications of the real samples hinder 
answering the initial question which motivated our 
research, namely whether the ground state of an isolated 
tetrahedron with antiferromagnetic exchange remains 
degenerate or whether this degeneracy is easily lifted by 
some small perturbation. Our study reveals that the 
behaviour of such " simple" systems is quite complex and 
is governed by high sensitivity of electronic structure of 
the system on the precise molecular configuration and 
this determines the magnetic behaviour of the material. 
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FIG. 15: (Color online) Vibrational density of states for each 
atom type (top). The vertical, dashed line highlights pro- 
nounced minima at 5 meV in all partial DOS (except for the 
oxygen vDOS) and the existence of lower frequency modes. 
Examples of latter for deuterium, carbon and copper, are 
shown separately (down). 
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FIG. 16: Two-atom (bond) and four-atom (twist angle) cor- 
relations as a function of frequency from the MD calculation. 




